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Abstract 

We discuss the phase diagram of the Shastry-Sutherland model for arbitrary 
spin S and derive rigorous lower and upper bounds on the phase boundaries 
of the dimer phase by using various versions of a variational ansatz in combi- 
nation with the exact diagonalisation method. 



I. INTRODUCTION 

We study a two-dimensional Heisenberg model with additional frustrating interactions 
of strength Ji on every second diagonal bond (see Fig. 1) given by the Hamiltonian 

H = J,Y^S-3; + J, SiSj, (1) 

<i,j> <i,j>Dimer 

where Si denotes a spin operator for spin S at site i. The model, commonly referred to as the 
Shastry-Sutherland model [1], displays a rich zero-temperature phase diagram as a function 
of Ji and J2, showing long range ferromagnetic, antiferromagnetic and helical order, as well 
as short range spin liquid behaviour. Also for large frustrating coupling Ji the model has 
an exactly known ground state, built up of uncoupled dimers on the diagonal bonds. It was 
the existence of this exact ground state, which has first attracted attention to the model. 
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Recently the Hamiltonian (1) for 5" = | has been widely discussed, because it represents 
essential features of SrCu2(B03)2, a newly discovered spin gap substance [2]. SrCu2(B03)2 
appears to have a ground state of dimer singlets and some of its unusual magnetic properties 
are well explained by the Shastry-Sutherland Hamiltonian. 



V 


1 








Y 























FIG. 1.: Shastry-Sutherland lattice with coupling J\ along the diagonal and J 2 along the 
vertical and horizontal bonds. 

In the classical limit {S — 00) of the Hamiltonian (1) the nature of the ground state 
phases and their boundaries are easily analyzed (see Fig. 2). The model shows two long 
range helical phases for < | J2I < Ji, which are separated by an antiferromagnetic dimer 
phase along the hne = J2 < Ji. In the regime IJ2I > J\ the ground state is ordered 
antiferromagnetically for J2 > and ferromagnetically for J2 < 0. For J2 = > Ji a phase 
of independent ferromagnetic dimers occurs. 

The ground state energy of the helical phase can be found by minimizing the energy of 
a cluster of only three spins. From this basic three spin entity the ground state of the whole 
system is constructed. Depending on the initial choice of spin directions one finds a spin 
helix with a four-fold degenerate direction and with a twist angle of ^ = tt ± arccos( J2/ Ji) 
between neighbouring spins [3]. 

Away from the classical limit quantum effects transmute the classical antiferromagnetic 
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FIG. 2.: Exact phase diagram of the classical Heisenberg model (left) and schematic phase 
diagram of the Heisenberg model for 5" > | (right). For S > ^ adjacent to the dimer phase 
two phases occur, whose nature is still a point of controversy. 

dimers for Ji > into unique singlet dimers and stabilize the resulting dimer phase. In what 
follows we will often use the variable x := J2/ Ji, which is the inverse of the frustration. For 
S < 00 and J2 > the dimer phase occupies a finite region of the phase diagram and at a 
certain critical value x"{S) > there is a phase transition of first order [3] to a new phase, 
whose nature is controversially discussed at the moment. For even larger x the system 
changes to the antiferromagnetic regime. 

In the ferromagnetic region (J2 < 0) the situation is slightly simpler, since ior S — ^ 
there is a first order phase transition directly from the dimer to the ferromagnetic phase 
at X — —1. This transition marks one of the exactly known points in the phase diagram. 
For S > ^, like in the antiferromagnetic regime, the system crosses at a negative xl{S) from 
the dimer to an intervening phase. Finally at x = —1 the ferromagnetic region is reached 
for all S. 
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For Ji < the classical phase diagram is not changed by quantum effects. The fer- 
romagnetic dimer phase for J2 — becomes a phase of independent spin-25' dimers with 
macroscopic degeneracy. It may be interesting to note, that the spin-^ square lattice an- 
tiferromagnet for = Ji < J2 is continuously connected to the spin-25' square lattice 
antiferromagnet as J2 > > Ji — > — 00. 

There are two major problems connected with the phase transition between the short 
range spin liquid phase and the adjacent phases. Firstly the nature of the phases adjacent 
to the singlet dimer phase - both on the ferromagnetic and on the antiferromagnetic side - 
is not clear and a point of intense investigation at the moment. Secondly for arbitrary S the 
boundary of the dimer phase cannot be determined by a simple argument like in the case 

In this paper we concentrate on the second problem and give rigorous lower (Sec. II) 
and upper (Sec. Ill) bounds on the phase boundaries of the dimer phase by using various 
versions of a variational ansatz for finite clusters. 

Various approximative methods have been used recently to determine this phase bound- 
ary. In Ref. [3] the phase boundary was calculated by means of Schwinger boson mean-field 
theory. By the fiow equation method in Ref. [4] a value of x'^{S = ^) = 0.63 was obtained 
and in Ref. [5] an intervening plaquette phase was suggested for = | occupying the regime 
0.861 > X > 0.677. In Ref. [6] by an extension of SU{2) to the symplectic groups Sp{2N) 
an intervening phase with helical and incommensurate order was found between a dimer 
phase and a region with collinear commensurate order. Also in Ref. [7] it was found for 
a two-dimensional model with frustration, that the dimer and the antiferromagnetic phase 
are separated by an intervening regime, which is characterized as a weakly incommensurate 
spin density wave. 
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II. LOWER BOUNDS ON x« AND ON -xl 



To obtain lower bounds on we decompose the Hamiltonian H of Eq. (1) into cluster 
terms with 

H = J2Hr (2) 

i 

and calculate the lowest eigenvalues E'}^ of the Hamiltonian Hj^ of the clusters (N is the 
number of spins in the cluster). The decomposition of H is chosen in such a way, that the 
clusters cover the whole lattice without overlapping bonds. Taking the ground state of H as 
variational state for the finite clusters [8] it follows that E^^ is always smaller than or equal 
to the ground state energy E^ of H. To obtain a bound on we choose finite clusters (see 
Fig. 3) which have the dimer state as an eigenstate and calculate their critical value x'^ j^. 
It is obvious, that x^^y is always smaller than the a;^ of the infinite system. 

As the simplest possible system we consider a plaquette with four spins (first entry in 
Fig. 3). In this case the energy of the second lowest state, which for growing x is competing 
with the dimer state, crosses the energy of the dimer ground state at 



From what was said above a;" 4 (5") is a strict lower bound for x". (Note that for 5" > 1/2 the 
first estimate Eq. (3) is better than the bound J1/J2 > 2(1 + S) given in Ref. [3].) Better 
limits on were obtained by calculating the ground state energies of the systems shown in 
Fig. 3 using the Lanczos method. 

For a system with 31 sites we thus obtain a best lower bound of 0.5914 for Xc{S — |). We 
want to point out, that even for 5 = |, where systems of up to 31 spins were calculated, it is 
difficult to make a good finite size analysis with the results of the finite clusters, because not 
only the system size, but also the shape of the cluster influences the ground state energy. On 
these grounds the significance of the recent estimate of Xc{S = |) = 0.7 ± 0.01 in Ref. [9], 
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which was obtained by an extrapolation of three systems with different stiapes, appears 
questionable. 

To show the tendency of a possible extrapolation we plot in Fig. 4 the results for Xc for 
5" = I as a function of 1/y/N and connect by a straight line the results of systems with 
N — 12, 24 spins which have corresponding shapes. The point Xc — 0.65, where this line 
meets the ordinate shows a reasonable value for Xc, but obviously cannot be taken as a 
serious extrapolation. The best bounds we obtained for 5" = 1, |, 2 are given in Fig. 3. 
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FIG. 3.: for 5* = |, 1, |, 2 calculated on finite clusters. The negative numbers denote x{ 
for >S = 1,|,2. 

Let us now consider the ferromagnetic regime. We find, that the dimer state of a four 
spin plaquette is the ground state for a; > — ^ which implies 
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4is) < (4) 

This bound is exact ior S — since x^. — —1 coincides with the boundary of the ferromag- 
netic phase. For 5" > | we again find an improvement on the result of the four spin system 
by considering larger clusters. Some results for 5* = 1, |, 2 are shown in Fig. 3. 
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FIG. 4.: x'^{S = I) for the systems shown in Fig. 3 plotted versus 1/\/N. The dotted line 
connects the results for N — 12, 24 sites, but should not be taken as an exact lower bound. 




III. UPPER BOUNDS ON AND ON -xl 

Upper bounds for the stabihty of the dimer phase can be obtained by considering a 
variational ansatz for finite clusters. We only briefly sketch the idea here, since it is widely 
used in the hterature (see e.g. Ref. [10] and references therein). 

The Hamiltonian is split into clusters without common spins and external bonds con- 
necting the clusters, 

-f^ = ^cluster + -f^&ond, (5) 

cluster 
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as indicated for four site clusters in Fig. 5. We use as variational state Hduster l^ci), where 
\^ci) is the ground state of -f^^y^^r- total spin of the clusters is zero in the ground 

state the expectation value of the external bonds vanishes 



( II ^cl\Hbond\ n *c/>=0. (6) 

cluster cluster 

We have calculated the ground state energies of H^^^^^^ for clusters with up to N=32 
spins and obtained upper bounds on by comparing them with the energy of the dimer 
state. 




FIG. 5.: Four site clusters (shaded grey) surrounded by external bonds. The product state 
built up of the eigenstates of the cluster Hamiltonian is used as variational ansatz to obtain 
an upper bound on x". 

For a four spin plaquette and arbitrary spin S we thus find 

x« < ^ (7) 

For stripe configurations of the type 2 x y the resulting for 5" = | are shown as a 
function of one over system size in Fig. 6. Note that results for two different types of stripes 
(as indicated in the figure) are given and that both types extrapolate to almost the same 



point 0.7126(1) for — > oo. A similar extrapolation for stripes of width four (see. Fig. 7) 
gives an upper bound of = 0.6955(4) for 5' = |, which is better than the bound obtained 
for stripes of width two. Again three different shapes converge to about the same result as 
N approaches infinity. The lowest value in Fig. 7 was obtained ior a, N — 32 system yielding 
the exact bound < 0.7050. An extrapolation of stripes of width two for S=l and systems 
of up to 16 spins gives a best upper bound of x'^{S = 1) < 0.618. 
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FIG. 6.: Extrapolation of the upper bounds obtained from stripe like configurations. Linear 
extrapolation of the last two points yields 0.7126 and 0.7127. 

In the limit S ^ oo the upper bound (7) derived from the four spin system does not 
provide the correct behaviour expected for the classical Heisenberg model. We can however 
find another bound, which yields a better upper limit for 5" > 1 by using a helical product 
state as a variational state. 

The ground state energy for 5" = oo in the helical phase [3] is given by 

I I J2 

Ehel = ^ -J^ (8) 

where the angle between two neighbouring spins is 6* = tt ± arccos(J2/Ji) and the length of 
the spins is normalized to 1. 
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FIG. 7.: Same extrapolation as in Fig. 6 for stripes of width four. Linear extrapolation of 
the last two points yields 0.6959, 0.6955 and 0.6953 for the three different shapes. 

We consider as variational ansatz a system of (quantum) spins polarized along the direc- 
tions of the classical spins. Comparing the expectation value of this state with the energy 
of the dimer state we find 




-(^'y + j;)s'<-JiS{s + i) (9) 

from which 

results as a criterion for the instabihty of the dimer state. Thus for x > 1/ the dimer 
state is no longer the lowest state, because the hehcal product state has a lower energy. 
Since the ground state energy of the helical state does not depend on the sign of J2 Eq. (10) 
also gives an upper bound on —x^ in the ferromagnetic regime. 




IV. SUMMARY 

In this paper we have discussed the phase diagram of the Shastry-Sutherland model for 
arbitrary spin S. In the regime Ji < | J2I the phase diagram for the quantum models {S < 00) 
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fully agrees with that of the classical model {S — oo). Quantum effects do, however, modify 
the classical phase diagram in the regime < | J2I < Ji- As emphasized already by Shastry 
and Sutherland [1] the dimer phase is strongly stabilized by quantum effects and exists in 
a regime xl{S) < J2/J1 < x'^{S) which grows with decreasing spin S. The nature of the 
intervening phases originating from the classical helical phases is not totally clear, but their 
regime shrinks as the dimer phase expands. For J2 < the intervening phase between the 
dimer and the ferromagnetic phase is certain to exist for 5" > 1/2, but completely disappears 
for S = 1/2. For J2 > exact statements on the existence of an intervening phase between 
the dimer and the Neel phase are not available, but in Ref. [5] rather strong evidence in 
favor of its existence for 5" = 1/2 is provided. 




FIG. 8.: Summary of bounds for the dimer phase in the Shastry-Sutherland model. 

The emphasis of the present work was put on the derivation of exact upper and lower 
bounds on the boundaries xl{S) and x'^{S) of the dimer phase. Our results on this topic are 
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summarized in Fig. 8. The dark grey area shows the region where the dimer state is certain 
to be the ground state and the white area is the region where the dimer state is certainly 
not the ground state. From the results presented in this paper it follows, that the exact 
boundary must be located in the light grey area. The black lines in Fig. 8 represent curves 
discussed in the text. The upper bounds for x'^{S) were obtained by a hehcal variational 
ansatz {^^): by a finite cluster ansatz (^^) and for 5* < 1 by extrapolating results for 
series of finite clusters (see Figs. 6 and 7). The lower bounds on x'^{S) were derived from 
four spin plaquettes (2;^:j::i) from the finite systems shown in Fig. 3 for 5" < 2. In the 
ferromagnetic regime we have a lower bound {-^) from the helical ansatz, which is shown 
together with the upper bounds obtained from the largest tractable finite clusters (Fig. 3). 
The intervening phase on the ferromagnetic side extends down to x — —1 (bottom of Fig. 8) 
where the ferromagnetic regime begins for all S. 
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